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We present the results of a precision computation of Bk with Wilson fermions. Simulations are performed at 
different lattice spacings, enabling continuum limit extrapolations. Two different twisted mass QCD (tmQCD) 
regularisations are considered for the computation of bare matrix elements. In both cases the relevant four-fermion 
operator renormalises multiplicatively. In one regularisation it is possible to perform the computation directly 
at the physical kaon mass value, thus avoiding extrapolations in the mass. Nonperturbative renormalisation is 
carried out using available Schrodinger Functional results. 



1. tmQCD for Bk 

The kaon mixing parameter Bk is currently 
the dominant source of uncertainty in unitarity 
triangle analyses of CP violation. It is therefore 
of utmost importance to have as precise a lattice 
QCD determination of this quantity as possible. 
Our project aims at a computation of Bk with 
Wilson fermions that brings all the systematics 
(apart from quenching) under control, in order 
to obtain a final result with an uncertainty of at 
most a few percent. 

The main sources of uncertainty in existing 
quenched computations with Wilson fermions are 
the limits to simulations imposed by the presence 
of exceptional configurations and the mixing un- 
der renormalisation of the parity-even part of the 
relevant four-fermion operator 
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Both problems can be eliminated with the use 
of tmQCD Wilson regularisations.^ A framework 

*TaIk presented by C. Pena at Lattice 2004, Fermilab. 
^For a review of tmQCD and recent developments related 
to it, see [3. 



to compute Bk without the need of determin- 
ing mixing coefhcients, via a tmQCD regularisa- 
tion with a fully twisted (w, d) doublet and an 
untwisted s quark, was introduced in [2]. The 
first results obtained within this approach have 
been presented in previous Lattice conferences [2]. 
Here we present our first analysis of the full set 
of results, and supply a preliminary value for Bk 
in the continuum limit. 

We always work within a Schrodinger Func- 
tional (SF) framework. Our action is nonper- 
turbatively 0(a) improved in the bulk of the SF 
cylinder, and one- loop 0(a) improved at the time 
boundaries. For the nonperturbative renormal- 
isation of the four-fermion operator we use the 
results in 

2. Results 

Simulations have been performed at three val- 
ues of the lattice coupling /3 and a number of 
different quark masses in the limit of unbroken 
SU(3)y symmetry Mg — M^, where Mi are the 
physical rcnormalised quark masses. The param- 
eters employed in the runs from which our final 
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value for is derived are summarised in Ta- 
bleQ] together with the values of the pseudoscalar 
meson mass and the bare value of Bj^ obtained in 
each case. The latter is extracted from a suitable 
ratio of correlation functions, as described in 0. 

To obtain the value of Bk at the physical point 
we have extrapolated the results in Table to 
(roMps)^ = 1.591, assuming a linear dependence 
of Bk on A/pg (which actually leads to our best 
fits). We stress that an independent simulation 
has been performed at each value of the mass, 
and therefore the results are fully uncorrelated. 

After nonperturbative renormalisation, we can 
extrapolate our results to the continuum limit as 
depicted in Fig.^ Since the matrix element is not 
fully 0(a) improved, the renormalisation group 
invariant parameter Bk is expected to approach 
the continuum limit linearly in a/rp; however, we 
observe no cutoff dependence within the statisti- 
cal accuracy of the results, and therefore take a fit 
to a constant as our best result. As a preliminary 
value in the continuum limit we finally quote 
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Bk = 0.817(22) , 

(2 GeV) = 0.592(16) . 

We have performed a number of tests to ensure 
that the quoted error adequately takes into ac- 
count a number of systematic uncertainties. First 
of all, simulations on larger lattices (correspond- 
ing to L « 2.2 fm) at /3 = 6.0 have shown that 
the value of Bk at the lightest pseudoscalar me- 
son mass does not depend on the physical vol- 
ume within its statistical uncertainty. We also 
performed a study, working at /? = 6.0 and 
roMpg ~ 1.777, of the dependence of Bk on the 
breaking of the SU(3)v symmetry, parameterised 
by e = {Ms - Md)/{Ms + Md). Fig.Hshows the 
results: no significant dependence is observed up 
to e ^ 0.4. 

3. Eliminating mass extrapolations 

The main remaining source of systematic un- 
certainty in the previous method comes from the 
fact that using an untwisted s quark prevents us 
from computing at physical values of the kaon 
mass in the Mg — Md limit. There are vari- 
ous tmQCD regularisations that allow to avoid 




Figure 1. Extrapolation (solid points, band) to 
the continuum limit of the renormalisation group 
invariant _B-parameter Bk- To allow for an easy 
comparison with a reference computation employ- 
ing Wilson fermions, we include the results of 
the method without subtractions of |S] (empty 
points), which display a linear extrapolation in 
a/ro. Also displayed (triangles) are the results 
for Bk obtained with (7r/4)-twisted quarks (see 
text). 
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Figure 2. Dependence of the bare Bk aX (3 — 6.0 
on the SU(3)y breaking parameter e. 



this problem. One simple possibility consists in 
considering a regularisation in which the doublet 
(s,d) is twisted with angle a — 7r/4. It can be 
easily shown that the renormalisation properties 
of the relevant four-fermion operator are exactly 
the same as in the previous case. This approach 
is designed for mass degenerate strange and down 
quark masses, and, since both quarks are twisted, 
the physical kaon mass can now be reached with- 
out any problem. An alternative that uses only 
fully twisted (a = 7r/2) quarks has been proposed 
in 6 . An attractive feature of this latter proposal 
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6.0 


16^ X 48 


402 


0.1335 


0.135169 


0.03186 


0.3895(11) 


1.025(8) 






398 


0.1338 


0.135178 


0.03152 


0.3551(12) 


0.993(11) 






402 


0.1340 


0.135183 


0.02708 


0.3320(11) 


0.977(11) 






400 


0.1342 


0.135187 


0.02261 


0.3054(11) 


0.943(10) 


6.2 


X 64 


200 


0.1346 


0.135780 


0.028324 


0.2815(8) 


0.980(7) 






201 


0.1347 


0.135783 


0.025985 


0.2685(9) 


0.952(8) 






214 


0.1349 


0.135787 


0.021290 


0.2432(10) 


0.934(8) 


6.3 


243 X 72 


204 


0.1348 


0.135771 


0.023639 


0.2385(10) 


0.970(14) 






212 


0.1349 


0.135773 


0.021254 


0.2251(11) 


0.937(12) 






205 


0.1351 


0.135776 


0.016467 


0.1964(12) 


0.885(15) 



Table 1 

Parameters for the runs from which our main data have been extracted. The bare twisted d quark mass 
has been denoted by fid- The mass parameters are tuned so as to have equal renormalised s and d quark 
masses up to O(a^) cutoff effects. 
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Figure 3. Interpolation to obtain Bk{P — 6.0) 
at the physical kaon mass in the (7r/4)-twisting 
approach (solid point) from the results for four 
different pseudoscalar masses (open points). 



is that it is expected to bring in automatic 0(a) 
improvement of the (bare) matrix element. 

As a feasibility study of the (7r/4)-twisting ap- 
proach we have performed simulations aX (3 — 
6.0,6.2 and a number of different pseudoscalar 
meson masses, on a box of physical length L « 
1.5 fm. The value of Bk at the physical kaon 
point is then obtained by linear interpolation in 
Mpg, as shown e.g. in Fig.|31 The results are also 
shown in Fig. ^ They exhibit remarkably small 
uncertainties, and suggest a good agreement with 
the outcome from the previous approach. How- 



ever, preliminary results with the same run pa- 
rameters on a larger physical volume hint at the 
presence of a small but, given the small statis- 
tical uncertainties, still relevant finite volume ef- 
fect. This issue has therefore to be settled by 
more extensive simulations. 
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